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ABSTRACT 

We report Fermi-LAT observations of the radio-loud AGN SBS 0846+513 (z=0.5835), 
optically classified as a Narrow-Line Seyfert 1 galaxy, together with new and archival 
radio-to-X-ray data. The source was not active at 7-ray energies during the first two 
years of Fermi operation. A significant increase in activity was observed during 2010 
October-2011 August. In particular a strong 7-ray flare was observed in 2011 June 
reaching an isotropic 7-ray luminosity (0.1-300 GeV) of l.OxlO 48 erg s _1 , comparable 
to that of the brightest flat spectrum radio quasars, and showing spectral evolution in 
7 rays. An apparent superluminal velocity of (8.2±1.5)c in the jet was inferred from 
2011-2012 VLBA images, suggesting the presence of a highly relativistic jet. 
Both the power released by this object during the flaring activity and the apparent 
superluminal velocity are strong indications of the presence of a relativistic jet as pow- 
erful as those of blazars. In addition, variability and spectral properties in radio and 
7-ray bands indicate blazar-like behaviour, suggesting that, except for some distinct 
optical characteristics, SBS 0846+513 could be considered as a young blazar at the 
low end of the blazar 's black hole mass distribution. 

Key words: galaxies: active - galaxies: nuclei - galaxies: Seyfert - galaxies: individual 
(SBS 0846+513) - gamma rays 



1 INTRODUCTION 

Active Galactic Nuclei (AGN) are the most luminous persis- 
tent sources of high-energy radiation in the Universe. How- 
ever, only a small percentage of AGN are radio-loud, and 
this characteristic is commonly ascribed to the presence of 
relativistic jets, roughly perpendicular to the accretion disk. 
Accretion of gas on to the supermassive black hole (SMBH) 
is thought to power these collimated jets, even if the nature 
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of the coupling between the accretion disc and the jet is 
still among the o utstanding open questions in high-energy 
astrophysics (e.g., iBlandford Il2000l ; Weier 1120031 ). Certainly 
relativistic jets are the most extreme example of the power 
that can be generated by a SMBH in the centre of an AGN, 
with app arent bolometric luminosities up to 1Q 49 ~ 5 erg 
s" 1 (e.g.. lAckermann et af]|2010l ; iBonnoli et alj|201lf ). and 
a large fraction of the power emitted in 7 rays. 

Before the launch of the Fermi satellite only two 
classes of AGN were known to generate these structures 
and thus to emit up to the 7-ray energy range: blazars 
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and radio galaxies, both hosted in giant elliptical galaxies 
l|Blandford fe Rees I ll97Sh . The first two years of observa- 
tions by the Large Area Telescope (LAT) on board Fermi 
confirmed that these two are the most numero us classes of 
identified sources in the extragalactic 7-ray sky (lAbdo et al.l 
2010a; Nolan et al.| [2012). but the discovery of variable 7-ray 
emission from 4 radio-loud Narrow-Line Seyfert 1 galaxies 
(NLSls) revealed the presence of a possible emergin g third 
class of AGN with relativistic jets (lAbdo et a l, 2009a, bjc]). In 
addition, Fermi-LAT observations gave us the opportunity 
to study in more detail particular subclasses of the estab- 
lished ty pes of 7-ray emitting AGN (e.g. Broad L ine Radio 
Galaxies, lAbdo et al.ll2010bl ; iKataoka et"ai1l201ll ). On the 
contrary, no ra dio-quiet Seyfert galaxie s were detected in 7 
rays until now (|Ackermann et aT]|2012h . 

NLS1 are a cl ass of AGN identified 
by lOsterbrock fe Pogge I (| 19851 ) and characterized by 
the following optical properties: narrow permitted lines 
(FWHM (H/3) < 2000 km s" 1 ) emitted from the broad line 
regio n, [OIIIJ/H/3 < 3 (a criterion introduce d by Goodrich! 
1 19891 ). and a bump due to Fe II (see e.g. I Pogge I I2OO0T 
for a review). They also exhibit strong X-ra y variability, 
steep X-ray spectra (ph oton indices Fx > 2; iBoller et al.l 
1996; Grupe et al, l20ld ). and substantial soft X-ray excess 
( Boiler et al. 1 119961 ). These characteristics point to systems 
with smaller masses of the central black hole (10 6 -10 8 M , 
lYuan et al.l |2008| ) than in blazars and radio galaxies and 
higher accretion rates (close to or above the Eddington 
limit, lYuan et al l |200S| ). NLS1 are generally radio-quiet 
(R < 10, radio-loudness R being defined as ratio of rest- 
frame 1.4 GHz and 4400 A flux densities), w ith only a small 
fract ion classified as radio-loud (< 7%, iKomossa et al.l 
2006), even more sparse (~2.5%) are very radio-loud 
NLSls (R > 100), while generally 10%-15% of quasars 
are radio- loud and very radio- loud. In the past, several 
authors investigated the peculiarities of radio- loud NLSls 
with non-simultaneous radio to X-ray data, suggesting 
similarities with t he young stage of quasars or different 
types of blazars ( Komossa et all 120061 ; lYuan et all 120081 : 
iFoschini et "all [2009). The strong and variable radio emis- 
sion, and the flat radio spectrum together with variability 
studies suggested the presence in some radio- loud NLSls 
of a relativistic jet, confirmed in some objects by the 
Fermi-LAT detection in 7 rays (|Abdo et al.ll2009d ). This 
finding poses intriguing questions about the nature of these 
objects, the onset of production of relativistic jets, and the 
cosmological evolution of radio-loud AGN. The impact of 
peculiar characteristics of the central engines in radio-loud 
NLSls, which seem quite different of those of blazars and 
manifest in their peculiar optical characteristics, on the 
7-ray production mechanisms is currently under debate. 

In June 2011 high 7-ray activity from SBS 0846+513 
was obse rved by Fermi-LAT. Pre liminary results were re- 
ported in lDonato fc Perkinsl (|201 lh . That was an important 
confirmati on of the detec tion of a new 7-ray NLS1 after the 
claim bv lFoschinil (|201ll fl. 



1 The analysis presented in Foschini i201ll ) was performed with 
the P6_V3 IRFs, starting from the first Fermi-LAT catalogue 
(1FGL) as a reference for the background sources, and over the 
period August 2008-February 2011. 



First identified with a BL Lac object at redshift z = 1.86 
during a bright state bv lArp et all (| 19791 ), SBS 0846+513 is 
a high redshift object positioned very clos e to a low redshift 
galaxy at the end of a chain of five galaxies. lArp et all (|l979l ) 
noted also a small nebulosity close to SBS 0846+513 that 
could be a normal red galaxy observed by chance near the 
object due to projection effects. The Sloan Digit al Sky Sur- 
vey (S PSS) spectrum of the source reported in IZhou et al.l 
(2005) is typical of NLS1: FWHM (H/3) = (1710 ± 184) 
km s" 1 , [OIII]/H/3 ~ 0.32, a nd a strong Fe II b ump. Sim- 
ilar values were obtained by lYuan et alj {2008) analysing 
the same SDSS spectrum: FWHM (H/3) = (1810 ± 191) km 
s _1 and [OIII]/H/3 ~ 0.31. Moreover, the SPSS sp e ctrum 
showed that its true redshift is z = 0.5835. iNottale 1 (|l986l ) 
and others suggested that it was a gravitationally lensed 
quasar with variability due to gravitational amplification by 
a star in an intervening galaxy. However, with broader wave- 
length coverage and high resolution, no signs of any inter- 
vening galaxies have been found in the SDSS spectrum. The 
Hubble Space Telescope (HST) image, collected when SBS 
0846+513 was in a faint state (V ~19.7 mag) does not show 
signific ant resolved struct ure, thus no indication of the host 
galaxy (|Maoz et al.lll993l ). Hig h polarization (> 10%) was 
det ected in optical and radio bv lMoore fc Stockman! (|l98ll) 
and lSitko et al.l l| 19841 ). A remarkable optical variability was 
observed in the past with an amplitude of AV ~ 5 mag ove r 
1 year and AV ~ 4 mag over 1 month l|Arp et al.lll979l ). 
The high pola rization, very hig h brightness temperature 
(T b > 10 13 K; IZhou et "all 120051 ) . and the large-amplitude 
variability observed in optical were clues for the presence of 
an at least mildly relativistic jet, now confirmed by the LAT 
detection in 7 rays and superluminal motion revealed by our 
new radio observations. 

In this paper we present the detection of SBS 0846+513 
in 7 rays and discuss its characteristics by means of the avail- 
able radio-to-7-ray data. Preliminary results about the high 
re solution VLBA obse rvations of the source were presented 
m lOrienti et all (|2012l ). The paper is organized as follows: 
in Section 2 we report the LAT data analysis and results. 
In Section 3 we report the result of the Swift observations 
performed in August-September 2011. Radio data collected 
by Medicina, OVRO, Effelsberg, VLBA, and VLA are pre- 
sented in Section 4 and discussed in Section 5 and 6. In 
Section 7 we present the spectral energy distribution (SED) 
modeling, while discussion and concluding remarks are pre- 
sented in Section 8. 

Throughout the paper, a A-CDM cosmology with Hq 
= 71 km s" 1 Mpc~\ fl A = 0.73, and fi m = 0.27 is adopted. 
The corresponding luminosity distance at z = 0.5835 is d_L = 
3406 Mpc, and 1 arcsec corresponds to a projected distance 
of 6.584 kpc. 



2 FERMI-LAT DATA: SELECTION AND 
ANALYSIS 

The Fermi-LAT is a 7-ray telescope operating from 20 MeV 
to > 300 GeV. The instrument is an array of 4 x 4 identical 
towers, each one consisting of a tracker (where the photons 
are pair-converted) and a calorimeter (where the energies of 
the pair-converted photons are measured). The entire instru- 
ment is covered with an anticoincidence detector to reject 
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Figure 1. Monthly integrated flux (E > 100 MeV) light curve of 
SBS 0846+513 obtained from 2010 August 4 to 2011 August 4. 
Arrows refer to 2-<r upper limits on the source flux. Upper limits 
are computed when TS < 10. The horizontal line indicates the 
period considered for the SED. 



the charged-particle background. The LAT has a large peak 
effective area (~ 8000cm 2 for 1 GeV photons), an energy 
resolution typically ~10%, and a field of view of about 2.4 
sr with an angular resolution (68% containment angle) bet- 
ter than 1° for energi es above 1 GeV . Furt her details about 
the LAT are given bv lAtwood et all (|2009h . 

The LAT data reported in this paper were collected over 
the first 40 months of Fermi operation, from 2008 August 4 
(MJD 54682) to 2011 December 4 (MJD 55899). During this 
time the LAT instrument operated almost entirely in survey 
mode. The analysis was performed with the ScienceTools 
software package version v9r23pl. The LAT data were ex- 
tracted within a 15° Region of Interest (Rol) centred at 
the radio location of SBS 0846+513. Only events belong- 
ing to the "Source" class were used. In addition, a cut on 
the zenith angle (< 100°) was also applied to reduce con- 
tamination from the Earth limb 7 rays, which are produced 
by cosmic rays interacting with the upper atmosphere. The 
spectral analysis (from which we derived spectral fits and 
photon fluxes) were performed with the post-launch instru- 
ment response functions (IRFs) P7S0URCE_V6 using an un- 
binned maximum likelihood method implemented in the Sci- 
ence tool gtlike. 

The background model used to extract the 7- 
ray signal includes a Galactic diffuse emission compo- 
nent and an isotropic component. The model that we 
adopted for the Galactic component is given by the file 
gal_2yearp7v6_v0.fits, and the isotropic component, which is 
the sum of the extragalactic diffuse emission and the resid- 
ual charged particle background, is parametrized by the file 
iso_p7v6source.tx10. The normalizations of both components 



in the background model were allowed to vary freely during 
the spectral point fitting. 

We examine the significance of the 7-ray signal from 
the sources by means of the Test Statistics (TS) based 
on the likelihood ratio test. The Test Statistic TS = 
2Alog(likelihood) between models with and without the 
source is a measure of the probability of having a 7-ray 
source at the localization specified, which compares models 
whose parameters have been adjus ted to maximise the likeli- 
hood of the data given the model jMattox et al.lll996t ). The 
source model used in gtlike includes all the point sources 
from the second Fermi-LAT catalogue (2FGL; iNolan et al.l 
I2012D that fall within 20° from the source. The spectra of 
these sources were parametrized by power-law functions, ex- 
cept for 2FGL J0920. 9+4441, and 2FGL J0957.7+5522, for 
which we used a log-parabola in their spectral modeling as 
in the 2FGL catalogue. We removed from the model the 
sources having test statistic TS < 25 and/or fluxes below 
l.OxlO" 8 photons cm" 1 s" 1 over 40 months and repeated 
the fit. Thus a final fitting procedure has been performed 
with the sources within 10° from SBS 0846+513 included 
with the normalization factors and the photon indices left 
as free parameters. For the sources located between 10° and 
15° we kept the normalization and the photon index fixed 
to the values obtained in the previous fitting procedure. The 
Rol model includes also sources falling between 15° and 20° 
from the target source, which can contribute to the total 
counts observed in the Rol due to the energy depended size 
of the point spread function of the instrument. For these ad- 
ditional sources, normalizations and indices were also fixed 
to the values of the 2FGL catalog. 

SBS 0846+513 was not in the 1FGL or 2FGL cata- 
logues, indicating that the source was not detected with TS 
> 25 in either one year or two years of Fermi observations 
l|Abdo et al.ll2010al ; INolan et al.ll2012T ). Integrating over the 
first two years of Fermi operation the fit yielded a TS = 
14. The 2-cr upper limit over the first 2 years of LAT data 
is 8.5xl0" 9 photons cm" 2 s" 1 in the 0.1-300 GeV energy 
rang«[f|, assuming the photon index is T = 2.3, which is the 
average value found for low-synchrotron- peaked blazars in 
the S econd LAT AGN Catalog (2LAC) (|Ackermann et all 
2012). On the contrary, the fit with a power-law model to 
the data integrated over the third year of Fermi operation 
(2010 August 4 - 2011 August 4; MJD 55412-55777) in the 
0.1-300 GeV energy range results in a TS = 653, with an 
integrated average flux of (6.7 ± 0.5) xlO" 8 photons cm" 2 
s" 1 and a photon index F = 2.23 ± 0.05. 

Finally, over the period 2011 August 4-December 4 
(MJD 55777-55899) only a TS = 8 was obtained, with a 2-cr 
upper limit of 1.9xl0~ 8 photons cm" 2 s" 1 in the 0.1-300 
GeV energy range, assuming the photon index is F = 2.23. 
During August-December 2011 the source was not detected 
with TS > 10 using a time bin of 1 month, indicating a 
significant variability in the 7-ray activity of the source on 
month time-scales. For the rest of the paper we focused on 
the 7-ray data collected over 2010 August 4-2011 August 4. 

The 7-ray point source localization by means of the 



3 To obtain a convergence we performed the fit considering only 
2 http://fermi.gsfc.nasa.gov/ssc/data/access/lat/Background the sources within 10° from SBS 0846+513 and E > 300 MeV 

Models.html and extrapolated the flux down to 100 MeV. 
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gtf indsrc tool applied to the 7 rays extracted during the 
third year of observation results in R.A. = 132.48 deg, Dec. 
= +51.13 deg (J2000), with a 95% error circle radius of 
0.06°, at an angular separation of 0.01° from the radio posi- 
tion of SBS 0846+513 (R.A. = 132.492 deg, Dec. = +51.141 
deg, J2000). This implies a strict spatial association with 
the radio coordinates of the NLS1 SBS 0846+513. The pro- 
posed association could be questioned to some extent due to 
the presence of a low redshift star-forming galaxy almost on 
the line of sight toward SBS 0846+513. And indeed, several 
members of the Local Group and also some nearby starburst 
galaxies have been recently detected by Fermi-LAT and thus 
established as sources of high energy photons. Their 7-ray 
luminosities range however between 10 40 and ~3 xlO 43 erg 
s _1 , and no variability in their 7-ray emis sion was neither 
observed nor expected l|Abdo et al.ll2010eh . The very large 
7-ray luminosity of the LAT object discussed here and its 
flaring activity therefore rule out the foreground galaxy as 
the source of the detected 7 rays. 

In order to test for curvature in the 7-ray spectrum 
of SBS 0846+513 two alternative spectral models with re- 
spect to the power-law (PL) were used: the log parabola 
(LP), dN/dE oc EIE~ a - p log ( B / B ^ (lLandau et all fl98r1 : 
iMassaro et alj 12004) . and the broken power law (BPL). In 
the case of the log parabola the parameter a is the spectral 
slope at the energy Eo and the parameter j3 measures the 
curvature around the peak. We fixed the reference energy 
Eo to 300 MeV. We used a likelihood ratio test to check the 
PL model (null hypothesis) against the LP model (alterna- 
tive hypothesis). T hese values may be compared, following 
iNolan et all l|2012t ). by defining the curvature Test Statis- 
tic TS C urve=2(log LLpTog Lpl) = 10 corresponding to ~3.3-<r 
difference. This value i s below the threshold of TS C urvc=16 
applied in INolan et all (|2012h for defining a significant cur- 
vature, thus suggesting only a hint of spectral curvature in 
the average 7-ray spectrum of SBS 0846+513. We tested also 
the broken power-law model, with photon indices Pi below 
and P2 above the break energy -EWcak- We fixed -Ebrcak at 1.4 
GeV. This value was estimated studying the profile of the 
Likelihood function, fixing -Ebrcak at different values between 
100 MeV and 3 GeV with a step of 50 MeV and finding the 
minimum log(likelihood) value. Also for this spectral model 
we obtain TS cur vc=2(log LepL-log Lpl)=10, therefore we 
adopt the power-law for the following LAT analysis. The fit 
results are reported in Table [1] 

Figure [T] shows the 7-ray light curve of the third year 
of Fermi observations built using 1-month time bins. For 
each time bin the photon index of SBS 0846+513 and all 
sources within 15° from it was frozen to the value resulting 
from the likelihood analysis over the entire year, while for 
the sources between 15° and 20° from SBS 0846+513 the 
photon index was frozen to the value reported in the 2FGL 
catalogue. If TS < 10 the value of the fluxes were replaced 
by the 2-cr upper limits. The systematic uncertainty in the 
flux is energy dependent: it amounts to 10% at 100 MeV, 
decreasing to 5% at 560 MeV, and increasing to 10% above 
10 GeV. By means of the gtsrcprob tool we estimated that 
the highest energy photon emitted by SBS 0846+513 (with 
probability >80% of being associated with the source) was 
observed on 2 May 2011 at distance of 0.34° from the source 
with an energy of 31 GeV. 

The source was not continuously detected over the en- 
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Figure 2. Integrated flux (E > 100 MeV) light curve of 
SBS 0846+513 obtained from 2011 June 4 to 2011 July 4 with 
3-day or 1-day time bins. Arrows refer to 2-cr upper limits on the 
source flux. Upper limits are computed when TS < 10. 



tire third year and the flux remained below 10~ ph cm" 
s , except in the period 2011 June 4-July 4 when an in- 
crease in flux by a factor of ~ 4 was observed. Considering 
the high activity of the source we extracted a spectrum over 
that period, obtaining a photon index F = 1.98 ± 0.05 and 
a flux of (24.4 ± 2.1) xlO" 8 photons cm" 2 s" 1 . A simi- 
lar spectral evolution has been already observed in other 
flat spectrum radio quasars (FSRQs) during a bright state 
l|Abdo et al.ll2010dh . A light curve focused on the period 
2011 June 4-July 4, extracted with the photon index fixed to 
T = 1.98, was built with 3-day and 1-day time bins (Fig.0. 
We used a 1-day time bin for the period with higher statis- 
tics. The peak of the emission was observed between June 
30 19:43 UT and July 1 19:43 UT, with a flux of (87 ± 16) 
xl0~ 8 photons cm -2 s _1 in the 0.1-300 GeV energy range. 
The peak 7-ray flux is a factor of ~13 higher with respect 
to the average flux estimated over the entire third year of 
observation. The 7-ray flare is characterized by a sharp in- 
crease in flux of a factor of ~2.5 in 1 day and return to the 
previous flux level in 1 day. 

For taking into consideration the possible influence of the 
choice of the To (the start time of the bin) used for build- 
ing the light curve on the determination of the 7-ray peak 
we shifted the To of 12 hours back and forth and recalcu- 
lated the daily flux, obtaining (53 ± 15) xlO -8 and (66 ± 
15) xlO -8 photons cm~ 2 s -1 , respectively. Finally, during 
the month of high activity, replacing the power-law with a 
log parabola or a broken power-law we obtain a TS CUT .„ e 
= 49, indicating a significant curvature. We noted that the 
curvature parameter j3 in the log-parabola spectral model 
increased as the flux increased (see Table [TJ. 



3 SWIFT DATA: ANALYSIS AND RESULTS 



The Swift satellite (|Gehrels et al.l [20041 ) performed two ob- 
servations of SBS 0846+513 on 2011 August 30 and Septem- 
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Table 1. Unbinncd likelihood spectral fit results. 





PL 


LP 


BPL 


Time Period (MJD) 


r 


a /3 


Ti r 2 


55412-55777 


2.23 ± 0.05 


1.93 ± 0.12 0.13 ± 0.05 


2.00 ± 0.09 2.70 ± 0.18 


55717-55747 


1.98 ± 0.05 


1.23 ± 0.19 0.30 ± 0.07 


1.51 ± 0.12 2.81 ± 0.21 



Table 2. Log and fitting results of Sim/4/XRT observations of SBS 0846+513 using a power-law model with Nu fixed to Galactic 
absorption. "Observed flux. 



Observation 


Net Exposure Time 


Photon index 


Flux 0.3-10 keV a 


Date 


sec 


r 


xl0~ 13 erg cm -2 s _1 


2011-08-03 


6738 


1.51 ±0.24 


9.0 ± 2.2 


2011-09-15 


5427 


1.43 ±0.28 


7.7 ± 2.9 



ber 15 for 6.8 ksec and 5.5 ksec, respectively. The observa- 
tions were performed with all three onboard intruments: the 
X-ray Telescope fXRT: lBurrows et all (|2005m . 0.2-10.0 keV) 
the U ltraviolet Optical Telescope (UVOT; iRoming et al.l 
ll2005l). 170-600 nm) a nd the Burst Alert Telescope (BAT; 
iBarthelmv et~ai1 (|2005h . 15-150 keV). 

The hard X-ray flux of this source is below the sensi- 
tivity of the BAT instrument for such short exposures and 
therefore the data from this instrument will not be used. 
Moreover, the source was not present in the Swift BAT 
58-month hard X-ray catalogue (jBaumgartner et al. 20101) 
and t he 54-month Palermo BAT catalogue ( Cusumano et al.1 
l2Q10h . 

The XRT data were processed with standard procedures 
(xrtpipeline vO.12.6), filtering, and screening criteria by 
using the Heasoft package (v6.11). The data were collected 
in photon counting mode in both of the observations, and 
only XRT event grades 0-12 were selected. The source count 
rate was low (< 0.5 counts s _1 ); thus pile-up correction was 
not required. Source events were extracted from a circular 
region with a radius of 20 pixels (1 pixel ~ 2.36"), while 
background events were extracted from a circular region 
with radius of 50 pixels away from the source region. An- 
cillary response files were generated with xrtmkarf , and ac- 
count for different extraction regions, vignetting and PSF 
corrections. We used the spectral redistribution matrices 
v013 in the Calibration database maintained by HEASARC. 

Considering the low number of photons collected (< 
200 counts) the spectra are rebinned wit h a minimu m of 
1 count per bin and the Cash statistics l|Cash I Il979l ) are 
used. We fit the spectrum with an absorbed power-law us- 
ing t he photoelectric absorption model tbabs (|Wilms et al.l 
2000), with a neutral hydr ogen column fixed to its Galac- 
tic value (2.91 xlO 20 cm" 2 ; iKalberla et al.1 l2005h . The low 
photon statistics prevents us from fitting the X-ray data 
with a more complic ated model than a simple power-law. 
We noted that also in iGrupe et all fcoicft . in which a large 
sample of NLSls observed by Swift/XRT was studied, the 
absorption column density was fixed to the Galactic value 
and the authors found that in the majority of NLSls the 
spectrum can be fit sufficiently well by an absorbed single 
power-law model with negligible intrinsic absorption. The 
fit results are reported in Table [2] The relatively hard X-ray 



spect rum with respect to the other NLSls fe.g. lGrupe et all 
|2010F ) could be due to the contribution of inverse Compton 
radiation from a relativistic jet, similarly to flat spectrum ra- 
dio quasars. A photon index 1.5-1.8 was observed in X-rays 
also for PMN J 0948+0022, the firs t NLS1 detected in 7 rays 
by Fermi-h AT l|Abdo et al.ll2009bl ; iFoschini et al.ll201ll l. 

In the past SBS 0846±513 was detected in X-rays only 
by ROSAT, with a 0.1-2.4 keV flux o f 2.7xl0~ 13 erg c m" 2 
s" 1 and a photon index Y = 1.77l ;^ l|Yuan et alj|2008h . As 
a comparison the flux observed on August 30 by Swift/XBT 
in the same energy range with F fixed to 1.77 is 3.7 ± 0.6 
x 10~ 13 erg cm~ 2 s~\ ~40% higher than the ROSAT obser- 
vations, adding evidence of the variability character of the 
source. 

In 2011 August the UVOT instrument took 5 frames in 
the wl band only, while four sequences with all v, b, u, wl, 
m2, and w2 filters were acquired in 2011 September. These 
data were processed with v6.10 of the Heasoft package and 
the CALDB release dated 2011 August 12. Source counts 
were extracted from a circular region of 5" radius centred 
on the source, while background counts were derived from a 
circular region of 10" radius in the source neighbourhood. By 
fitting the source spectrum of September 15 with a power- 
law, we calculated the effective wavelengths, count rate to 
flux conversion factors (CFa), and Galactic extinctions for 
the UVOT bands accord ing to the procedure explained in 
iRaiteri et ail |20ld . l201ll '). The results are shown in Tabled] 
To obtain the UVOT de-reddened fluxes reported in Table [4] 
we multiplied the count rates for the CFa and corrected for 
the corresponding Galactic extinction values Aa- 

We note that a Swift/WOT observation performed on 
2011 September 15 found SBS 0846+513 about 0.4 magni- 
tude brighter in V-band with respect to the HST observation 
in January 1992. 



4 RADIO DATA: ANALYSIS AND RESULTS 

New and archival radio data for SBS 0846+513 were col- 
lected from Medicina, OVRO, the Effelsberg 100-m, VLA 
and VLB A. In the following subsections we present the ob- 
servations performed by these facilities. 
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Table 3. Results of the UVOT calibration procedure: effec- 
tive wavelengths A e ff , count rate to flux co nversion factors CFa , 
and Galactic extinction calculated from the lCardelli et ah I jl98Sh 
laws. 



Filter 


^eff 


CF A 


A A 




A 


10~ 16 erg cm -2 s -1 A -1 


mag 


V 


5430 


2.60 


0.09 


b 


4360 


1.47 


0.11 


u 


3476 


1.65 


0.14 


uvwl 


2621 


4.40 


0.19 


uvm2 


2257 


8.36 


0.24 


uvw2 


2087 


5.97 


0.23 



Table 4. Results of the analysis of UVOT data for SBS 0846+513. 



MJD 


Filter 


Magnitude 


Magnitude_error 


55803.644 


uvwl 


19.57 


0.07 


55819.497 


V 


19.10 


0.30 


55819.489 


b 


19.67 


0.19 


55819.488 


u 


19.22 


0.17 


55819.497 


uvwl 


19.41 


0.15 


55819.501 


uvm2 


19.55 


0.15 


55819.493 


uvw2 


19.55 


0.10 



4.1 Medicina 

We observed SBS 0846+513 with the Medicina radio tele- 
scope five times between 2011 August and November. We 
used the new Enhanced Single-dish Control System (ESCS) 
acquisition system, which provides enhanced sensitivity and 
supports observations with the cross scan technique. All ob- 
servations were performed at both 5 and 8.4 GHz; the typical 
on source time is 1.5 minutes and the flux density was cali- 
brated with respect to 3C 286. Since the signal to noise ratio 
in each scan across the source was low (typically ~ 3), we 
performed a stacking analysis of the scans, which allowed 
us to significantly improve the signal to noise ratio and the 
accuracy of the measurement. We list the final values of the 
5 and 8.4 GHz flux density in Table [5] 

4.2 OVRO 

As part of an ongoing blazar monitoring program, the Owens 
Valley Radio Observatory (OVRO) 40 m radio telescope has 
observe d SBS 0846+513 at 1 5 GHz regularly since the end 
of 2007 l|Richards et alj|201ll ). This monitoring program in- 
cludes over 1500 known and likely 7-ray loud blazars above 

Table 5. Results of the Medicina 32-m radio observations. 



Obs. date 


s 5 


Ss.4 




(mjy) 


(mjy) 


2011-08-10 


170 ± 15 


235 ± 20 


2011-09-08 


210 + 20 


185 ± 20 


2011-09-22 


190 ± 20 


180 ± 20 


2011-10-13 


190 ± 20 


280 ± 30 


2011-11-16 


250 ± 30 


275 ± 25 



2009 2010 2011 2012 




56000 



line [MJD] 



Figure 3. 15 GHz radio light curve for the period 2008 August 
4-2011 December 12 from the OVRO telescope. The downward 
arrow indicates the time of the peak of the 7-ray activity observed 
by Fermi-LAT. 

declination —20°. The sources in this program are observed 
in total intensity twice per week with a 4 mjy (minimum) 
and 3% (typical) uncertainty. Observations are performed 
with a dual-beam (each 2.5 arcmin FWHM) Dicke-switched 
system using cold sky in the off-source beam as the refer- 
ence. Additionally, the source is switched between beams 
to reduce atmospheric variations. The absolute flux density 
scale is calibrated using obser vations of 3C 286. adopting 
the flux density (3.44 Jy) from lBaars et all l| 19771 ). This re- 
sults in about a 5% absolute scale uncertainty, which is not 
reflected in the plotted errors. 

4.3 Effelsberg 100-m 

The centimeter spectrum of SBS 0846+513 was observed 
with the Effelsberg 100-m telescope on 2011 April 30 (MJD 
55681.8) within the framework of a Fermi-related mon- 
itoring program of 7 -ray blazars (F-GAMMA program; 
iFuhrmann et al.ll2007h . The measurements were conducted 
with the secondary focus heterodyne receivers at 2.64, 8.35, 
14.60, and 32.00 GHz. The observations were performed 
quasi-simultaneously with cross-scans, that is, slewing over 
the source position, in azimuth and elevation directions, with 
adaptive numbers of sub-scans for reaching the desired sensi- 
tivity (for details, see lFuhrmann et al . 2008; [Angelakis et al.l 
2008). Pointing offset correction, gain correction, atmo- 
spheric opacity correction, and sensitivity correction have 
been applied to the data. 

4.4 VLA and VLBA data 

To implement the information on the flux density variability 
available from the single-dish observations and to study the 
source structure on kpc and pc scales, we analysed archival 
VLA and VLBA data at different frequencies obtained be- 
tween 1986 and 2011. To study the proper motion we also 
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Table 6. Log of the archival VLA observations and flux density. 



Freq 


Date 


Code 


Obs.time 


Beam 


Flux 


GHz 






min 


arcsec X arcsec 


mjy 


1.4 


1986-04-10 


AV127 


9.16 


1.40x1.12 


194±6 


4.8 


1986-04-10 


AV127 


9.16 


0.44x0.35 


286±9 


4.8 


1996-01-05 


BW021 


1.66 


1.39x1.19 


332±10 


4.8 


1996-12-30 


BA018 


3.33 


0.43x0.36 


363±11 


4.8 


2009-05-09 


CALSUR 


0.66 


2.12x1.18 


196±6 


8.4 


1995-08-13 


AM484 


0.5 


0.24x0.22 


356±11 


8.4 


1995-09-02 


AM484 


0.5 


0.32x0.23 


327±10 



made use of 15-GHz VLBA data from the MOJAVE pro- 
graming. 

Usually SBS 0846+513 was observed with the VLA as a 
phase calibrator and the on-source observing time is gener- 
ally quite short. Logs of VLA and VLBA observations are 
reported in Tables [6] and [7] respectively. 
The data reduction of both VLA and VLBA data was per- 
formed following the standard procedures implemented in 
the National Radio Astronomy Observatory (NRAO) AIPS 
package. In the case of VLA data, the accuracy of the am- 
plitude calibration was checked by means of either 3C48 
or 3C 286, and it resulted to be within 3%. For the VLBA 
datasets, a priori amplitude calibration was derived using 
measurements of the system temperature and the antenna 
gains. Uncertainties of the amplitude calibrations were esti- 
mated to be within 5-10%. For MOJAVE data we imported 
the fully calibrated uv datasets (|Lister et al.ll2009h . 
The final VLA and VLBA images were produced after a 
number of phase-only self-calibration iterations. The source 
parameters, like flux density and angular size, were derived 
on the radio image plane by means of the AIPS task JMFIT 
which performs a Gaussian fit to the source and its subcom- 
ponents (Table [HJ. For the extended component visible at 15 
GHz, we determine the flux density by means of TVSTAT, 
which performs an aperture integration on a selected region 
on the image plane. 



5 RADIO VARIABILITY 

Multi-epoch studies of the radio emission of SBS 0846+513 
show significant flux density variability. In Fig. [3] we report 
the OVRO observations obtained since the launch of Fermi. 
From the analysis of the OVRO light curve at 15 GHz it 
is clear that high and low activity states are interspersed. 
Interestingly, no 7-ray emission was detected during the ma- 
jor radio flaring episodes that occurred before 2011. A max- 
imum flux density variation of a factor of ~2.0 was observed 
over 3 months during the period 2010 August-2011 August, 
while in the period 2008 August-2010 August a variation of 
a factor of ~2.8 was observed over a similar time interval. 
Flux density variability is detected also at 4.8 GHz, where 
the flux density varies about 70% during the time inter- 
val spanned by the VLA data (1986 to 2009). However, the 
poor time coverage of the VLA data does not allow us to 

4 The MOJAVE data archive is maintained at 
http://www.physics.purdue.edu/MOJAVE 
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Figure 4. Radio data of SBS 0846+513 collected at 5.0 GHz 
(triangles) and 8.4 GHz (squares) by Medicina, and 15 GHz (open 
circles) by OVRO during the period 2011 July 28 - November 19 
(MJD 55770-55884). 



fully characterize the flux density behaviour on such a long 
time interval. 

The detection of the 7-ray flare in 2011 June, triggered a 
monitoring campaign with the Medicina telescope. In Fig. [4] 
we report the Medicina flux density at 5 (triangles) and 8.4 
GHz (squares) in addition to the OVRO 15-GHz data (open 
circles) collected between 2011 July 28 and November 19. 
From August 20 (MJD 55793.628), i.e. about one month af- 
ter the 7-ray flare, the flux density at 15 GHz started to 
increase reaching the maximum of about 280 mjy around 
October 18, and then it started to decrease. At lower fre- 
quencies the flux density variation is observed after a longer 
time, firstly at 8.4 GHz (October 13), and then at 5 GHz 
(November 16). This time delay may be explained by opacity 
effects, which are more severe at longer wavelengths. Dur- 
ing this time interval also the shape of the radio spectrum 
changes, as is clearly visible in Fig. [S] (left panel). We must 
note that Medicina and OVRO observations are not strictly 
simultaneous and to build the radio spectrum we considered 
the OVRO observation that is the closest in time to the 
Medicina observing epoch. For comparison, in Fig. [5] (right 
panel) we show the simultaneous radio spectrum between 
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Figure 5. Left panel: Radio spectra of SBS 0846+513. Observations at 5 and 8.4 GHz are from the Medicina single-dish telescope, those 
at 15 GHz are from OVRO. Triangles, squares, and crosses refer to observations performed around 2011 August 10, October 13, and 
November 16, respectively. Right panel: Radio spectrum of SBS 0846+513 obtained by the Effelsberg single-dish on 2011 April 30 from 
2.64 to 32 GHz. 



VLA 8.4 GHz 
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Figure 6. VLA image at 8.4 GHz of SBS 0846+513. On the im- 
age we provide the peak flux density, in mjy/beam, and the first 
contour intensity (f.c, in mjy/beam) that corresponds to three 
times the noise measured on the image plane. Contour levels in- 
crease by a factor of 2. The beam is plotted on the bottom left 
corner of the image. 

2.64 and 32 GHz obtained by Effelsberg on 2011 April 30, 
i.e. before the high activity state detected in 7 rays. 



6 RADIO MORPHOLOGY 

The resolution of the VLA is not adequate to resolve the 
radio structure of SBS 0846+513 (Fig. [6j even when the 
array is in its most extended configuration (full width 



half maximum ^0.3 arcsec). When imaged with the high 
spatial resolution of the VLBA the source is resolved in two 
components with a core-jet structure (Figs. [7] and Et, a s 
also pointed out in a previous work bv lTavlor et al.l ([2005) . 
The core (component W) is still unresolved with an upper 
limit of 0.3 mas, while the jet structure (component E) is 
3.5 x 1.6 mas. The flux density ratio between components 
E and W is ~19 and 14 at 5 and 8.4 GHz respectively. 
Their separation is about 3.5 mas, corresponding to ~23 pc 
projected distance, given the source's redshift. At 15 GHz 
component E shows an extended low-surface brightness 
structure with a steep spectrum. On the other hand 
component W is resolved into two compact components 
(labelled Wl and W2 in Fig. |9]). To investigate a possible 
proper motion of the jet, we compared the separation 
between Wl, considered the core region, and W2, assumed 
to be a knot in the jet, at the three observing epochs. To 
this purpose, in addition to the data analysis on the image 
plane, we model-fitted the visibility data using gaussian 
components of the three-epoch MOJAVE data by means 
of the model-fitting option in difmap. This approach is 
preferable in case we want to derive small variations in 
the source structure and little changes in the position of 
the source components. From this comparison we found 
that Wl and W2 are separating with an apparent velocity 
of (0.38+0.07) mas/yr (Fig. I10p . which corresponds to 
(8.2±1.5)c. This apparent superluminal velocity suggests 
the presence of boosting effect. However, the availability of 
only three observing epochs spanning a short time baseline 
without frequent time sampling implies large uncertainties 
on the estimated values. 

When we compare the total flux density derived from 
VLBA data with that from the VLA, we find that the 
VLBA can recover at most 85% of the VLA flux density, 
even when the observations were performed simultaneously 
(i.e. 1996 January 6). This may be related to the slightly 
different observing frequencies, i.e. 4.8 GHz at the VLA, 
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Figure 7. VLBA image at 5.0 GHz of SBS 0846+513. On the 
image we provide the peak flux density, in mJy/beam, and the 
first contour intensity (f.c., in mJy/beam) that corresponds to 
three times the noise measured on the image plane. Contour levels 
increase by a factor of 2. The beam is plotted on the bottom left 
corner of the image. 



Figure 8. VLBA image at 8.4 GHz {bottom) of SBS 0846+513. 
On the image we provide the peak flux density, in mJy/beam, and 
the first contour intensity (f.c, in mJy/beam) that corresponds 
to three times the noise measured on the image plane. Contour 
levels increase by a factor of 2. The beam is plotted on the bottom 
left corner of the image. 



Table 8. Flux density in mjy of the components of SBS 0846+513 
from VLBA data. Epochs of the observations are reported in Ta- 
ble [7] 





E 


W 


Wl 


W2 


S5GHz,epl 


14±2 


269+19 






S5GHz,ep2 


17±2 


265+18 






S8.4GHz 


14±2 


196+14 






Sl5.3GHz,epl 


3.4+0.4 


192+15 


140+10 


52+4 


Sl5.3GHz,ep2 


3.1+0.3 


240+24 


175+17 


65+6 


Sl5.3GHz,ep3 


3.6+0.4 


209+20 


151+15 


58+5 



whereas it is 5.0 GHz at the VLBA. However, in this case 
the spectral index would have to be extremely steep (a r ~4; 
f v tx v~ ar ) to obtain such a difference in the flux density, 
making this explanation unlikely. Another possibility is 
related to observational limitations due to the lack of short 
spacings of the VLBA array, implying that only structures 
smaller than ~40 and 30 mas at 5 and 8.4 GHz, can be 
detected. This suggests that the missing flux on the parsec 
scale may be related to extended, low-surface brightness 
features like a jet component resolved out by the VLBA 
array. 
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7 SED MODELING 

The 7-ray detection of the radio loud NLSls has given us 
the possibility of studying the characteristics of this class of 
objects by modeling their broad band spectra. 

We have built a non-simultaneous SED, all of which is 



Figure 9. VLBA image at 15.3 GHz of SBS 0846+513. Data 
are from the MOJAVE programme. On the image we provide the 
peak flux density, in mJy/beam, and the first contour intensity 
(f.c, in mJy/beam) that corresponds to three times the noise 
measured on the image plane. Contour levels increase by a factor 
of 2. The beam is plotted on the bottom left corner of the image 
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Table 7. Log of the VLBA observations analysed in this paper and total flux density. a : data from the MOJAVE programme. 



Freq 


Date 


Code 


Beam 


Obs time 


Flux 


GHz 






mas X mas 


min 


mjy 


5.0 


1995-03-15 


BW015 


3.99x1.34 


50 


284±20 


5.0 


1996-01-06 


BW021 


4.13x1.29 


50 


281±20 


8.4 


2011-03-12 


BC196J 


2.00x0.97 


10 


210±16 


15.3 


2011-05-26 


BL149DP 


0.77x0.59 


30 


244±17 


15.3 


2011-07-15 


BL149DM a 


0.76x0.60 


30 


213±15 


15.3 


2012-01-02 


BL178AE a 


0.74x0.64 


36 


231±21 



0.8 
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Figure 10. Changes in separation with time between components 
Wl and W2. The solid line represents the regression fit to the 15- 
GHz MOJAVE data. 



based on non-flaring data. This ought to adequately repre- 
sent the emission from this object in a low state. The LAT 
spectrum was built with data centred on 2010 October 4 to 
2011 June 4 (MJD 55473-55716). In addition we included 
in the SED the Effelsberg radio data collected on 2011 April 
30 and the Swift (UVOT and XRT) data collected on 2011 
September 15, thus well before and after the 7-ray flare. The 
data from 2MASS All-sky PSC source catalogue and WISE 
preliminary release Source Catalogue provided information 
about the IR part of the spectrum. The flare centred on ~ 
MJD 55750 had a variability timescale of ~ 1 day, which 
constrains the size of the emitting region during the flare. 

We modeled the SED assuming emission from a rel- 
ativistic jet with mechanisms of synchrotron, synchrotron 
self-Compton (SSC), and Compton-scattering of a dust 
torus external to the jet (EC-dust ) . The des cription of the 
model can be found in lFinke et all |2008T l and [5 ermer et al.l 
(2009). Jet powers were calculated assuming a two-sided jet. 
Although the flaring timescale does not constrain the emit- 
ting region during the low state, we have produced a model 
fit which is roughly consistent with this timescale (about 1 
day). The synchrotron component can adequately explain 
the IR through UV points, although some of the 2MASS 
points are significantly higher than the rest, and are not 
well-fit. This is probably due to contamination by IR emis- 
sion from the host galaxy. We note that this IR excess could 
be consistent with the typical starlight expected from an el- 
liptical galaxy. The radio spectrum is flat in F v (Fig. [5| and 
probably the result of a superposition of several jet compo- 



nents (e.g., lKonigi"lll981i ), so for our purposes these points 
are considered upper limits. The synchrotron component we 
use is self-absorbed below ~ 10 12 Hz. An attempt was made 
to fit the X-ray through 7 ray data with an SSC component, 
but this was not found to be possible. Instead, an external 
component was required. Correlations of 7-ray and optical 
flares with radio light curves and rotations of optical polar- 
ization angles in low-synchrotron-peaked blazars seem to in- 
dicate the 7-ray/optical emitting region is outside the broad 
line region, where the dust torus is the likely seed photon 
source (e.g., iMarscher et~alll2010r ). So we chose as our seed 
photon source a dust torus, which was simulated as a one- 
dimensional ring with radius Rdust aligned orthogonal to the 
jet, emitting as a blackbody with temperature Tdust and lu- 
minosity Ldust- Observational evidence for a dust torus in 
NLSls as well as in the othe r Seyfert galaxies was recently 
reported in iMor et al.l (|2012T ). The fit can be seen in Fig. 1111 
and the parameters can be found in Table [H A descr iption 
of the parameters can be found in lDermer et al.l <|2009h . The 
dust parameters were chosen so that Rdust is roughly con - 
sistent with the sublimation radius (|Nenkova et aljfeoOfit ). 
assuming the torus luminosity is about 1/10 of the disc lumi- 
nosity. If the emitting region takes up the entire width of the 
jet, and the jet is assumed to be conical, its half-opening an- 
gle would be about 1°, approximately consist ent with those 
meas ured from VLBI observations of blazars (jjorstad et al.l 
2005). This size scale is roughly the same as for the 15 GHz 
core (Fig. |9}, and so if this model fit is correct, the 7-ray 
emitting region is probably located near this region, which 
is presumably a synchrotron self-absorption photosphere. 

The electron distribution used, a broken power-law with 
index pi — 2.2 below the break at ^' brk and p2 = 3.2 for 
Ibrk < 7' j i s consistent with particles injected with index 
2.2, and emission taking place in the slow cooling regime 
(e.g., iBottcher fc Dermerl 2002). That is, particles are in- 
jected between "/' min and 'y'max, with a cooling break at 

3m e c 2 
Ibrk = - 

iCaTlescUtot 

where utot is the total energy density in the blob frame, 
which in the case of our model fit is dominated by the ex- 
ternal energy density. The escape time t esc m R' b /c so that 
for our model fit 7^ rfc ~ 300. Also note that in this model fit 
the magnetic field and electrons are nearly in eq uipartition. 

Using the SDSS Data Release (DR) 1, IZhou et al.1 
( 2005) estimated the black hole mass from several methods, 
the H/3 broad line, the [OIII]A5007 narrow line, and the host 
galaxy's bulge luminosity, as ~ 8.2 xlO 6 M , ~ 5.2xl0 7 M Q , 
and ~ 4.3 x 10 7 M©, respectively. Large uncertainties (0.4- 
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Figure 11. Spectral energy distribution data (circles and 
squares) and model fit (solid curve) of SBS 0846+5f3 with the 
model components shown as dashed curves. The data points were 
collected by Effelsberg (2011 April 30), Swift (UVOT and XRT; 
2011 September 15), and Fermi-LAT (2010 October 4-2011 June 
4), together with archival data from 2MASS and WISE. The LAT 
spectrum during the flaring period (2011 June 4— July 4) is shown 
as the diamonds. 



0.7 d e x) are asso c iated to these techniques (jVerstergaardl 
12004 ) . IShen et al.1 (1201 ll ) estimated the BH mass based on 
the R/3 broad line and Mgll line in the SDSS DR7 spectrum 
as log(A/ SH /M ) = 7.99 ± 0.12 and log(M SH /M Q ) = 7.79 
± 0.16, respectively. Note that the broadband SED shows no 
evidence for a blue bump, so the optical spectrum probably 
has considerabl e contamina t ion fr om jet emission not taken 
into account in Sh en et al.l (|201ll ). and s o these mass esti - 
mates should be taken with caution. If the lZhou et all (|2005l ) 
H/3 estimate is correct, the total jet power will exceed the 
Eddington luminosity for this source ( Lsdd = 1-3 x 10 45 
erg s" 1 M B h/(W 7 M©) ). Thus the BH mass is probably 
more on the high side of these estimates, if the jet power is 
Eddington-limited. 

The Compton dominance for SBS 0846+513, i.e., the 
ratio of the peak luminosities of the Compton and syn- 
chrotron components, is ~ 7, which is a rather standard 
value for FSRQs, although quite a high value for BL Lacs 
(e.g., iGiommi et ai1l201ll ; iFinke I [2012 '). Its SED in many 
ways resembles that of a FSRQ, with an X-ray spectral in- 
dex ~ 1.5, consis tent with values for FSRQs in the BAT 
|Aiello et al.ll2009l ) and BeppoSAX (|Donato et al.ll2005l ) cat- 
alogues. As we discuss in section|S] its LAT flux and spectral 
index are also consistent with values for FSRQs. 

The flaring LAT spectrum is also shown in Fig. 1111 Un- 
fortunately, there are no multiwavelength data simultaneous 
to this flare. For FSRQs, 7-ray flares are often associated 
with optical flares, although there are occasions when they 
are not. Both types of flares are seen, for ex ample, from the 
FSRQ PKS 1510-089 l|Marscher et al.ll2010l ) so we can only 
guess at the behavior of SBS 0846+513 during this flare. 
However, note that during the flare, the object reaches a 
luminosity of L lAso ~ 1.0 x 10 48 erg s _1 , making the radia- 
tive power L 7i ,. ad ~ L ltiso /T 2 ~ 4.4x 10 45 erg s _1 , assuming 
T = 15. This is about 1/2 of the Eddington luminosity for a 
10 7 M black hole. 



8 DISCUSSION AND CONCLUDING 
REMARKS 

After the 4 o bjects detected by F ermi-LAT in the first year 
of operation |Abdo et al.1 l2009d ), SBS 0846+513 is a new 
NLS1 detecte d by Fermi-LAT during high 7-ray activity 
in 2011 June (|Donato &: Perki"nsll2011 ). The power released 
by this object during the flaring activity was a strong in- 
dication of the presence of a relativistic jet as powerful as 
those in blazars, supported by the apparent superluminal 
velocity of the jet derived by tracking the position of a jet 
component in 2011-2012 VLBA data (see Sec. [6]). Before 
the 7-ray flaring episode, the simultaneous multifrequency 
observations performed by Effelsberg showed a flat radio 
spectrum up to 32 GHz. After the flare, the spectral shape 
changed, becoming convex. The spectral variability was also 
accompanied by variations in the radio flux density, which 
were originally detected at the higher frequencies, later mov- 
ing to lower frequencies, likely due to opacity effects. These 
spectral and variability properties indicate blazar-like be- 
haviour , which has already b een observed in other 7-ray 
NLSls ((Fuhrmann et al.ll201ll ). 

iGhisellini et al.1 (|2008l , l201ll ), investigating the 7-ray 
properties of blazars detected by Fermi-LAT in the first year 
of observation, suggested a transition between BL Lac ob- 
jects and FSRQs that can be justified mainly by the different 
accretion regimes: highly sub-Eddington in the former case, 
near-Eddington in the latter. In this context the high accre- 
tion rate of SBS 0846+513, and more generally of the 7-ray 
radio- loud NLSls, is another indication of its similarity with 
FSRQs. A comparison of the BL L acs and FSRQs fro m the 
First LAT AGN Catalog (1LAC; lAbdo et al.ll2010d ) with 
"misaligned AGN" detected by Fermi (MAGN; non-blazar 
AGN with jets pointed away from the observer) in the T 7 - 
L 7 plane has shown that MAGN and blazars occupy differ- 
ent regions of the plane, with only two high redshift FRII 
galaxies, 3C 207 an d 3C 380, which lie among the FSRQs 
|Abdo et ai1l2010bh . This is in agreement with the idea that 
MAGN are less beamed than blazars. In this context it is 
interesting to consider also the 7-ray NLSls. For a direct 
comparison with the results shown in lAbdo et all |2010bl ) 
we calculated the flux and photon index over the third year 
of Fermi operation in the 0.1-10 GeV energy band, resulting 
in T 7 = 2.19 ± 0.06 and Flux (0.1-10 GeV) = (6.6+0.6) 
x 10 -8 erg cm -2 s _1 . The corresponding observed isotropic 
7-ray luminosity in the 0.1-10 GeV is 3.6xl0 46 erg s _1 . We 
plotted these values of SBS 0846+513 in the F 7 - L 7 plane 
toge ther with the FSRQ s, BL Lacs, and misaligned AGN 
from lAbdo et al.l |2010bl ). As can be seen SBS 0846+513 lies 
in the blazar-region, in particular in the transition region be- 
tween the distribution of BL Lacs and FSRQs (Fig. [T2T) . We 
note that also the other 7-ray NLS1 observed in flaring activ- 
ity, PMN J0948+0022 (F 7 = 2.26 ± 0.08 and 0.1-10 GeV 
luminosity of 9.6 xl0 46 erg s -1 over the first 24-month of 
Fermi operation; lGrandill201ll ), occupies the same blazar- 
region in that plane. This should reflect a similar viewing 
angle with respect to the jet axis and beaming factors for 
the 7-ray emission between blazars and the two 7-ray NLSls 
SBS 0846+513 and PMN J0948+0022. In the same way the 
spectral evolution during the flaring activity in June 2011 
observed in 7 rays from SBS 0846+513 is a common be- 
haviour in bright FSRQs and low-synchrotron-peaked BL 
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Table 9. Model parameters for the SED shown in Fig, 1111 
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Lacs detected by Fermi (|Abdo et al 1 l2010dh . with a change 
in photon index < 0.2-0.3 and an increasing spectral curva- 
ture. 

One of the key questions is the maximum power re- 
leased by the jets of NLSls. During the 7-ray flaring ac- 
tivity observed in 2011 June-July SBS 0846+513 reached 
an observed isotropic 7-ray luminosity (0.1-300 GeV) of 
l.OxlO 48 erg s _1 on daily timescales, compar able to that of 
lumin ous FSRQs. After PMN J0948+0022 ijFoschini et all 
l201ll L this is the second NLS1 observed to generate such a 
high power. This could be an indication that all the radio- 
loud NLSls are able to host relativistic jets as powerful 
as those in blazars, despite the lower BH mass; alterna- 
tively some NLSls could have peculiar characteristics al- 
lowing the development of these relativistic jets. The mech- 
anism at work for producing a relativistic jet is not clear, and 
the physical parameters that drive jet formation is still un- 
der debate. One fundamental parameter could be the black 
hole mass, with only large masses allowing for the efficient 
formation of a rel ativistic jet. It was fo r example noted 
by some authors (|McLure fc Jarvisl |2004| ; iLiu et all 120061 ; 
ISikora et al"1l2007M that quasars with M BH > 10 8 M Q reach 
radio loudness 3 orders of magnitude greater than quasars 
with A/bh < 3 x 10 7 M Q . The large radio loudness of SBS 
0846+513 co uld challenge this id ea if the black hole masses 
estimated by IZhou et al.l |2005h are confirme d. According 
to the "modified spin paradigm" discussed in ISikora et al.l 
|2007l L another fundamental parameter for the efficiency of 
a relativistic jet production should be the BH spin, with 
SMBHs in elliptical galaxies having on average much larger 
spins than SMBHs in spiral galaxies. This is due to the fact 
that spiral galaxies are characterized by multiple accretion 
events with random orientation of angular momentum vec- 
tors and small increments of mass, while elliptical galaxies 
underwent at least one major merger with large matter ac- 
cretion triggering an efficient spin up of the SMBH. The 
accretion rate (thus the mass) and the spin of the black hole 
may be related to the host galaxy, leading to the hypothe- 
sis t hat relativistic je ts can develop only in elliptical galaxy 
(e.g. lMarscher1l2009l ). 

In this context the discovery of relativistic jets in a 



class of AGN usually hosted by spiral galaxies, such as 
the NLSls, was a great surprise. Unfortunately only very 
sparse observations of the host galaxy of radio- loud NLSls 
are available at this time, and the sample of objects studied 
by iDeo et all ^OOtj ) and IZhou et all (|2006l ) had z < 0.03 
and z < 0.1, respectively, including both radio-quiet and 
radio-loud objects. We note that BH masses of radio-loud 
NLSls are generally larger with respe ct to the entire sampl e 
of NLSls (M BH » (2 - 10) x 10 7 M Q iKomossa et ailbOOrj ). 
even if still small when compared to radio-loud quasars. The 
larger BH masses of radio-loud NLSls with respect to radio- 
quiet NLSls could be related to prolonged accretion episodes 
that can spin-up the BHs. The small fraction of radio-loud 
NLSls with respect to radio- loud quasars could be an indi- 
cation that in the former the high accretion usually does not 
last s ufficiently long to significantly spin-up the BH (|Sikora I 
l2009h . 

Of the 7-ray emitting NLSls, host galaxy imaging is 
available only for 1H 0323+342, with HST and the Nordic 
Optical Telescope. These observations reveal a one-armed 
galaxy morphology or a circumnuclear ring, sug gesting two 
possi bilities: the spiral arms of the host galaxy (Zho u et al. 
1 20071) or the residual of a merging galaxy (|Anton et al. 
2008). No significant resolve d structures have been observed 
by HST for SBS 0846+513 jMaoz et al.lll993l ). and no high- 
resolution observations are available for the remaining 7-ray 
NLSls. Thus the possibility that the development of rela- 
tivistic jets in these objects could be due to strong merger 
activity, is not ruled out. Further high-resolution observa- 
tions of the host galaxies of 7-ray NLSls will be fundamental 
to obtain important insights into relativistic jet formation 
and development. 

To conclude, SBS 0846+513 shows all the characteris- 
tics of the blazar phenomenon. The extreme power released 
by SBS 0846+513 during the high 7-ray activity in 2011 
July confirms that, as with PMN J0948+0022, NLSls can 
host relativistic jets as powerful as blazars. Radio and 7- 
ray data collected for SBS 0846+513 suggest spectral and 
variability properties similar to blazars, and the modeling of 
the average SED gives similar results to those of blazars, in- 
cluding similar Lorentz factors. This could be an indication 
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Figure 12. The photon index T 7 of FRI radio galaxies (red cir- 
cles), FRII radio galaxies (green squares), BL Lacs (open blue cir- 
cles) and FSRQs (open black squares) are plotted together with 
the NLS1 SBS 0846+513 (magenta point including a black trian- 
gle) as a function of the obser ved isotropic 7-ray l uminosity (100 
MeV-10 GeV). Adapted from lAbdo et al.1 l|2010bl l. 
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that these 7-ray NLSls are low mass (and possibly younger) 
version of blazars. The detection of new radio-loud NLSls 
in 7 rays by Fermi-LAT will be important for extending the 
source sample, and better characterizing this new class of 
7-ray emitting AGN. Equally important will be to perform 
further multifrequency observations of the 7-ray emitting 
NLSls already detected by Fermi and investigate in detail 
their characteristics over the entire electromagnetic spec- 
trum, to help understand their nature. 
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